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Abstract 

We compute the form factors parametrizing radiative leptonic decays of heavy 
mesons B + — > r ye + i / for photon energies much larger than Aqcd, where per- 
turbative QCD methods for exclusive processes can be combined with the 
heavy quark effective theory. The form factors can be reliably obtained in 
this region in an expansion in powers of A/E^. The leading term in this ex- 
pansion displays an additional spin symmetry manifested in the equality of 
form factors of vector and axial currents. The leading twist form factors can 
be written as the convolution of the B meson light-cone wave function with 
a hard scattering amplitude, which is explicitly calculated to one-loop order. 
The Sudakov double logarithms of the form (— log — j^-) n are resummed to all 
orders. As an application we present a method for determining the CKM ma- 
trix element | V u b | from a comparison of photon spectra in B and D radiative 
leptonic decays. 
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Typeset using REVTgX 



I. INTRODUCTION 



The radiative leptonic decay B + — ► has received a great deal of attention in the 

literature as a means of probing aspects of the strong and weak interactions of a heavy 

quark system. The presence of the additional photon in the final state can compensate for 
the helicity suppression of the rate present in the purely leptonic mode. As a result, the 
branching ratio for the radiative leptonic mode can be as large as 10~ 6 for the /i + case ||, 
which would open up a possibility for directly measuring the decay constant /b |. A study 
of this decay can offer also useful information about the CKM matrix element \V u b\. 

Preliminary data from the CLEO collaboration indicate an upper limit on the branching 
ratio B(B + -> ^e + v) of 2.0 x KT 4 at the 90% confidence level jfo]] . With better statistics 



expected from the upcoming B factories, the observation and experimental study of this 
decay could become soon feasible. It is therefore of some interest to have a good theoretical 
control over the theoretical uncertainties affecting the relevant matrix elements. 

The hadronic matrix element responsible for this decay can be parametrized in terms of 
two formfactors defined as 

= <7(P7>e)lMl-75)?|S(u)>= (1) 



E(p,£*,v,p y )fv(Ey) + ^[e*(w -p 7 ) - (p 7 ) M (e* ■ v)]f A (E 7 ) . 

The photon energy in the rest frame of the B meson is E 7 = v-p 1 . The absolute normalization 
of the matrix element ([!]) can be fixed, in the limit of a soft photon, with the help of heavy 
hadron chiral perturbation theory. In the limit of a massless final lepton, which we will 
consider everywhere in the following, the leading contributions to fy{E 1 ) come from pole 
diagrams with a J p = 1~ intermediate state, and those to /a(£/ 7 ) from J p = 1 + states [|TJ. 
The dominant contribution comes from the B* + state, which is degenerate with B + in the 
heavy mass limit 

where the ellipses stand for contributions from higher states with the same quantum numbers 
and A = ins* — rriB- Qq is the light quark electric charge in units of the electron charge. 
The hadronic parameter j3 ~ 3 GeV -1 parametrizes the BB*^ coupling in the heavy quark 
limit ||. The formfactors fv,A have been also computed in the constituent quark model 
@H, using light-cone QCD sum rules || and in a light-front model f7j. 

Momentum conservation in the B + —>■ 7z^£ + decay can be written as rngt) = q + p-y, with 
q being the momentum of the lepton pair. The photon energy is given by 

E ^ = V -^ = — '2^- B (3) 

and depending on the invariant mass of the lepton pair, < q 2 < m 2 B it takes values within 
the window < £ 7 < In this paper we study the radiative leptonic decay B + — > ^ue + 

in the kinematical region A QCD where the perturbative QCD methods developed in 



TT| for exclusive processes can be applied. 
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In the limit m& — > oo and the photon energy satisfying E 1 3> Aqcd, the large momentum 
of heavy quark is carried away by the lepton pair and does not affect the hadronic part of 
the decay. Therefore, it becomes convenient to subtract away the large m^v component of 
the b quark momentum and define a new transfer momentum q by 

q = q + m b v (4) 

and g M = C(m°). In the kinematical region E 7 ^> Aqcd this momentum is space-like 
q 2 = A 2 — 2KE 1 < with A = tub — m& being the binding energy. 

Introducing the subtracted momentum q one notices that in the leading rrib — > oo limit 
the kinematics of our problem is very similar to the one for the tt° — > 7(jo 7 ) + J*(q) decay 
discussed in ]TI[. This suggests to apply QCD factorization theorems in order to expand 



the form factors fy^Ey) in inverse powers of q 2 , or equivalently 1/E 1 . In this formalism 
the form factors of interest can be written as the convolution of a hard scattering amplitude 
with the transverse momentum dependent wave function of the B meson, i[>(k+, k±). 

We will work in a reference frame where the photon moves along the "— " light-cone 
direction and has light-cone components of the momentum p 7 = (0, 2E 7 /v + , OjJ{] The 
transfer momentum is given by q = (Av + ,A/v + — 2E 7 /v + , 0j_). Then, to leading order in 
1/Ey (leading twist) and a s we find 

urn - me,) = q^i. J + o m . (5, 

The wave function ip(k + ,k±) depends on the "+" light-cone component, k + , and transverse 
momentum, k±, of the light quark momentum in the B meson. Its properties are studied in 
Sec. 2, where its moments are related to matrix elements of local heavy-light operators. The 
expression ([5]) for the form factors fy^{E^) as integrals over the light-cone wave function 
i/j(k + ) is derived in Sec. 3. The radiative corrections to this result induce a logarithmic 
dependence on E J} in addition to the power law \jE 1 . These include doubly logarithmic 
Sudakov corrections and mass-singular logarithms of the light quark mass log(m), which 
are resummed in Sec. 4. A few numerical estimates made with the help of a model wave 
function are presented in Sec. 5, where we present also a method for extracting the CKM 
matrix element \V u b\ from a comparison of the photon spectra in B and D radiative leptonic 
decays. A few details concerning the calculation of the radiative corrections are presented 
in an Appendix. 



II. LIGHT-CONE B MESON WAVE FUNCTION 

We consider a heavy B meson with the flavor content bq and momentum p = tjibV moving 
along the z axis. Its light-cone wave function can be expanded into a sum of multiparticle 
Fock components \B) = \bq) + \bqg) + • ■ -. The valence component is written explicitly as 



Throughout the paper we shall use the following definition of the light-cone components fc„ = 
(k + , k-, kj_) with k± = ko ± k^ and kj_ = (k±, fo)- 
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\B) = £ k ± )-^(\q a (k, ])b\k\ {)) - |g a & [)b\k\ ]))) . (6) 

The light quark g and heavy quark b in the B meson have light-cone momenta k = (k + , k±) 
and rribV + k', respectively. This gives the constraints k + + k' + = Av + and k± + k' ± = 0, 
with A = ins — rrib the binding energy of the B meson. The range of variation of k + is the 
interval (0, (m^ + A)v + ), corresponding to k' + = (Av + , —rribV + ). The wave function ip(k + ) 

only takes values significantly different from zero for k + /v + ~ A. 

The light-cone wave function ip(k + , k±) is related to the usual Bethe-Salpeter wave func- 
tion ty a/3 at equal light-cone "time" r = x° + x 3 

= J d 4 ^ k <(0\T + hp(0)P(0,Oq*(O\B(m B v)) (7) 



as 



-> r°° dk 

■* a p(k+,k ± )= ^ a p(k) (8) 

J-oo An 

= i, ^ C ,, ==; ' ~m ( u »(k, t)i>/j(v, 1) - u a (k, i)vp(v, t)) ip{k + , k±) . 
sj{2m b v + ){2k + ) V2 V J 

The quark fields appearing in the definition of the Bethe-Salpeter wave function are quan- 
tized on the light-cone 

q a (x) = £ (a- kjX u a (k, AJe"** + b[ x v a (k, A)e**) , (9) 

fe,A 

where the creation and annihilation operators satisfy {a~ k A , a~, ,} = 5^ fcSx^r. The light-cone 
spinors are defined as in |12| and are normalized according to u(k, X) , y + u(k, A) = 2k + . 



The static heavy antiquark field h(x) is related to the usual b field by b(x) = e mbV ' x h(x) 

and satisfies i/h{x) = —h(x). The path-ordered factor P(0, £) = Pe dzA ^ i s introduced 
to ensure gauge invariance of the Bethe-Salpeter wave function. 



Using the explicit expressions for the light-cone spinors u(k, A') and v(k', A) given in \1 I 
one finds the following result for the wave function (|S[) in the limit of an infinitely heavy b 
quark 

9ap{k+, k±) = --^Li/,(k+, fcj.) { (K + a ± ■ fc ± )A + l±^ 75 ) , (10) 
V2v + k + { 2 ) ap 

satisfying the usual on-shell conditions 

jM(k+, k±) = 0, k ± W= fcj.) . (11) 

We denoted with A + = 7_7 + /4 the projector on the space of fast-moving particles along the 
+z axis. 
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It is convenient to define the one- dimensional wave function^ ip{k + ) by integrating over 
the transverse momenta 

^(k+) = j ^^(k+,k ± ) (12) 

which satisfies 

1 [n~ c f A 1 + 1/.. | , r de-_*^/a,„.™: ^„/„ /! 

a/3 

Multiplying both sides with (7+75)^ gives 



V > i^tM ^ = /°° ^e^-/ 2 (0|T^(0) J P(0,O^(^-n_)|i?(p)> . (13) 



2iV c V#+) = jT ^ e ^«-/ 2 (0|T/ i (0)7 + 75g(2e^-)| J B(m B t;)) . (14) 

This relation can be used to express the moments of ip{k + ) in terms of matrix elements 
of local operators. To see this, the time-ordered product on the right-hand side of ([[4]) is 
expanded into a power series of the separation on the light-cone. This gives 

y/2NMk+) = £ (^-) neifc+ ^ /2 ( l^(^ + ) n g|5(m^)) . (15) 

Taking the j th moment with respect to k + one obtains the desired connection to local 
operators (see also ||T3| ) 



JNJ2 | o °° dk + (k + )^(k + ) = (0\h +l5 (tD + yq\B(m B v)) . (16) 



The first few moments of the wave function can be simply expressed in terms of known 
hadronic quantities. For j = the corresponding matrix element on the RHS of (|16"D is 
determined by the decay constant of the B meson in the static limit defined as 

(0|/i7M75<?|£( m £^)) = fBrriBVf, . (17) 
One finds the normalization condition 

■.M£M£) = /J w (i8) 

The first moment j = 1 is given by the matrix element 

f d (£) * (y = ffwh + ^D +)q \B M ) = f gay™ ■ da) 



2 We denote both ip(k) and ip(k+) with the same letter. The distinction between them is made 
through their arguments. 
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This result agrees with the intuitive notion that the averaged spectator quark momentum is 
proportional to the binding energy of the heavy hadron. To prove it, one starts by writing 
the most general form for the following matrix element, compatible with Lorentz covariance 



{Q\h^^ b (iD v )q\B(m B v)) = ag^ + bv^v v 



(20) 



The equation of motion for the light quark field ipq{x) = implies the constraint 4a + 6 
Another equation for these parameters can be obtained with the help of the relation 



0. 



Av u (0\h ltMl5 q\B(m B v)) = (0\h^ 5 (iD v )q\B) + (0\h^ 5 (iD v )q\B) 



(21) 
0, 



Multiplying both sides with v v and using the static quark equation of motion iv ■ Dh(x) 
one obtains a + b = v^A/em^. Solving for a and b gives the result presented in (|i~9|). 

In the presence of radiative corrections, the connection between the light-cone wave 
function and matrix elements of local operators is changed. For example, the zeroth moment 
(|T8|) acquires a scale dependence typical of matrix elements of operators in the effective 
theory with heavy quarks [T^JT^]. In the rest frame of the B meson (v + = 1) this is given 
by, after MS renormalization, 




-rrfB{^h)m B 



(22) 



dk. 



d 2 k A 

W) 



:ip{k+, kj 



1 + 



Air 



3 m 

--log — + A IR (k )+F(k ) 
. 1 



The term Am(ko) contains an IR singularity, which is regulated with dimensional regular- 
ization in D = 4 — 2e dimensions. The quantities Ajfj(ko) and F(ko) are given by 



A IR (k 



F(k Q 




(23) 



(24) 



We denoted here N^ R = \ — + log(47r) and k = v ■ k = \{k + + |^). The IR singularity 
in A/#(fco) originates from soft-gluon exchange between b and u in the initial state. The 
coefficient of N^ R depends on the angle between the momenta of the b and u quarks 
cosh d = (v ■ k)/m and is well known as the QCD bremsstrahlung function. Notice that it 
receives an imaginary contribution due to the instantaneous (Coulomb) interaction. 



The scale-dependent parameter fsi^h) m (|22|) is related to the physical decay constant 

f B by ujrjj 
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*-fesr«*>- (25> 

The logarithmic dependence on fi^ on the right-hand side of ( P^ ) matches that of the pa- 
rameter fsif^h), as it should. The remaining mass-singular logarithm can be absorbed into 
the wave function by introducing a factorization scale A 2 satisfying m <C Aq C d C A < /i h . 
Writing log(m 2 / /i 2 ) = log(m 2 /A 2 ) + log^ 2 //!^), the first term is absorbed into the wave 
function and the second is resummed into a factor similar to the one in (p5|) P|- 

The IR singular terms can be resummed to all orders in a s using the QCD evolution 
equations. In the resulting expression the contribution of multiple virtual soft gluon emis- 
sions exponentiates and it can be factorized out from the wave function. This suggests to 
absorb the IR singular term Air (as well as the mass singular logarithm log(A 2 /m 2 ) as 
explained above) into the light-cone wave function. For the purpose of normalization alone, 
one can define thus a modified wave function to one-loop order 

$(k+, k±,A) = ij(k + , k ± ) exp ( ^A m (v ■ k) + 0(a 2 s )) ( 1 + ^ • ? log (26) 



., An ' K " s/ J \ ' An 2^°m 2 7 

satisfying the normalization condition 

~2~, fa s (m b )\ 2/Po „ fd 2 k ±7/7 , A ./ <x s Gp . ,\ , s 

N c hmB {ir4) = L ^/^''^rf*')' (27) 

It will be shown below that the hard scattering amplitude to one-loop order is IR finite only 
when convoluted with this modified wave function. 




III. LEADING TWIST ANALYSIS OF B -> 7^ 

To leading order in a s there are two diagrams contributing to the matrix element ([[]), 
shown in Fig. 1. Only the diagram (a), where the photon is emitted from the light line, 
contributes to leading order in l/m^. Using the wave function (M), it can be written as 

r " = «• / W Tr K' iM1 - (28) 

The trace can be easily computed with the help of the basic relations 

Tr {(&;+ + a± ■ k ± )A + P v >y 5 7,375} = -v + (k + n + + k±)p , (29) 
Tr{(A; + + a_L • kx)A+P v ^ 5 7/3} = -~v + ie(ri-, k±, n+, (3) , (30) 



3 Note that this redefinition of the wave function applies strictly for the purpose of the normal- 
ization condition. The precise A-dependence of the wave function ifi(k+,A) is derived below in 
Sec. IV. 
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to which the expression (p8|) can be reduced by application of the identity 



l^lalv = g^alv + Qvalu ~ 9^1 'a + V, 0t, 0)7/375 . (31) 

We find in this way the following results for the form factors to tree level 

The formfactor fo^E*/) of the tensor current is encountered when considering the radiative 
rare decay B — > z/z/7. It is defined by the matrix element 

-^==(-f( Pl ,e)\b(T^ 5 q\B(v)) = i/ T (£ 7 )[(p 7 ) M e* - e*(p 7 )„] . (33) 

(The tensor structure e*v v — i^e* is forbidden by gauge invariance.) The matrix element 
of the ba^q current can be obtained from this one with the help of the identity a^ v = 

The corrections to the result ([32]) arising from the coupling of the photon to the heavy 
quark (Fig. 1(b)) are suppressed by A/m&. In fact the leading term in this expansion is 
calculable in terms of known quantities only. The corresponding correction is given by 

{h) f dk + d 2 k ± ( t -m h j- g + j ' + m b \ 

Performing the trace gives for the heavy quark contribution to the form factors 

m) = Q b M-±=r r dk fj X ^ k ± )(l + o(±,^-))= Q b I^L (35 ) 



where we used the normalization condition (0) for the wave function. This correction is 
potentially important for the case of charmed meson decays. 

The equality of the form factors in ( ]3"2"D to lowest order in a s can be understood as the 
consequence of a larger symmetry group of the Green functions in Fig. 1 to leading order in 
1/E T To see this, one notes that the momentum of the light quark entering the weak vertex 
contains a large light-like component p = k — p 7 = —E y ri- + k, with n± = (1,0,0, ±1). 
Therefore a natural description of this quark is in terms of the light-cone component of the 
quark field q_ defined as 

q _ = e ^(n--x) A _ g ^ A = J+^t_ ; ( 36 ) 

satisfying A_g_ = g_ or fi_q = 0. The corresponding Dirac action reads, when expressed in 
terms of this component [[17] 



q(iP)q = g_(-zn_ • £>)?_ + 0(1/ , (37) 

which contains an additional SU(2) symmetry group compared with the original one. This 
can also be seen in terms of the Feynman rules for the light quark line 



S 



Pr ° pagat0r: i -2E,(J 6 ..k)l^ + ie = l -2(n^k) +i e + (38) 
vertex: -igij a = -ig(-n_)^t a + 0(1/ . (39) 

To leading order in 1/E^ and l/m&, the weak current bTq can be written as h^Tq_, with 
h$' the static b quark field satisfying 7°/i4 = —h^>. Using the properties of the fields g_ 
and one can derive the following relation 

h { vh»q- = -(n-)J$q- + ie 0o rt(n-) a h®'yprYsq- • (40) 

Taking the matrix element of ( flUD between (7(p 7 ,e)| and \B(v)), and noting that 

(>y( Pl ,e)\h®q\B(v)) = 0, gives 

(7(^,6)1^7^-15^)) = ^/3^-)«(7(^-,e)l^ ) 7/ 3 7 5 g-|5(t;)) , (41) 

which reduces to fy(E y ) = /a(S 7 ) in the rest frame of v. Note that this is very different 
from other symmetry groups appearing in particle physics like flavor or spin as it is not 
apparent in the hadron spectrum; rather it is a symmetry of an internal part of a Feynman 



diagram mediating a decay process. Similar arguments have been used in [[16] to derive 
relations among semileptonic form factors in B — > tc, p using the additional symmetry of the 
so-called large energy effective theory for the final state hadron \T7\ . 



In the following we will show by explicit calculation to one-loop order that the equality 
fv(E 7 ) = f^(E 7 ) is preserved beyond tree level, for the leading terms in an expansion of 
these form factors in powers of Radiative corrections change the simple power law 

1/E^ by introducing a logarithmic dependence on the photon energy. To leading order in 
1/E-y and l/m& these corrections are given by (with % = V,A,T) 

tfip\ f ^nb) Y 2/Po n [N~c 1 r AU t/j(k+) ( a s C F \ 
The first factor accounts for the different renormalization of the weak current in the static 



quark effective theory and QCD [jUJIS]. The dependence on the hybrid renormalization 
scale /ih cancels between this factor and the hard gluon correction in the effective theory i^. 

We consider in the following the one-loop radiative corrections to the diagram in Fig. 
1(a). The heavy-light vertex correction shown in Fig. 2(a) has the form (for a general weak 
current bTq) 

A, = ~ig 2 C F {[r] (E« - (v • g) JW) - [TQ- W-Tv ■ (q - k)]((l - x) JW>} . (43) 
The scalar integral is defined by (the definition of ((1 — x) J^) is given in the Appendix) 

j(«0 = f — (u) 

J {2ny {-v ■ I + v ■ k> + ie)((l + q - k) 2 + it){l 2 + it) ' 1 ; 

The heavy quark can be taken on-shell such that its residual momentum k satisfies v ■ k' — 0. 
In fact the integrals Jw and ((1 —x)j( a >) are free of infrared and collinear divergences. The 
exact results for these integrals are presented in the Appendix in (|A7|) , ( |A8|) . In the limit 
kj r /E 1 — > they have the asymptotic expansions 
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J(a) 



((l-x)J (a) ) 



(4vr) 2 £ 7 



(4vr) 2 ^ 



1 2 2£ 7 
- log - — - 

2 k-L 



3 



1 2 2-E 7 2£L 
-log -— +log-— 



(45) 
(46) 



The UV divergent integral £(°' is evaluated using dimensional regularization in Z) = 4 — 2e 
dimensions. One obtains 



E (a) 



(4tt) 



iV. 



UV 



log 



2£LA; 



/'i 



+ 2 



(47) 



with iVf v = 1/e - 7s + log(4vr). 

Combining these results one finds the following contributions from the heavy-light vertex 
correction to the li factors from the diagram in Fig. 2(a) 



t)(a) 
-V 



-A 



n? v - io K 2Kfh * 



log 2 —=L + log —=L - — 

n,<_l_ ru_L O 



+ 2. 



(48) 



The light-light vertex correction shown in Fig. 2(b) introduces a correction to the photon 
coupling of the form 



ft -> f + 2z# 2 CV { ^ [-(2 - 2e) j| 6) + 2(p 7 ■ fc)(jW - J [b) + J [b) - J< 



(49) 



+ 2(e*.k)^[-J^ + j[ b) -4 b) + 4 b) }} 



The scalar factors J^ 1 are defined by 

d 4 / 



(2tt) 4 [(Z +p 7 - k) 2 - m 2 + ie] [I 2 - 21 ■ k + ie] (I 2 + ie) 

^1] = "I" J A ^k/^ku + t/g ^ (kyP^ v + P^^kv) + Jg ^P-yfiP-yu ■ 



(50) 

(51) 
(52) 



These integrals have collinear singularities, which will be regulated by giving the light quark 
a mass m. Their explicit results in the limit m 2 <C p 7 • k are given in the Appendix (see 
Eqs. ( |A9|) ). The term proportional to (e* • k)yf in the vertex correction (f49|) vanishes after 

the integration over k±. Keeping only the first term amounts to a multiplicative correction 
of the lowest order result. Using the results Eq. ( |A9| ) one obtains the following contributions 
to the li coefficients from the diagram in Fig. 2(b) 



'V 



-A 



N uv _ i 2E^k + , n 2£ 7 A; + 



+ 2 log ■ 



- 1 



(53) 



The self-energy correction on the internal light quark line (Fig. 2(c)) contributes 

1 . (54) 
Finally, the box diagram (Fig. 2(d)) is given by 



?(<0 
■v 



A 



_ N uv + lQg 2R 1 k ± 
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B 



19 F J {2-Ky{-vl + ie){{l + k-p ry ) 2 -m 2 + ie){l 2 + 2l-k + ie){l 2 + ie) ' l ' Mj 



The term of order Z° in the loop momentum has an IR singularity, which is regulated as 
before using dimensional regularization. The total contribution of the box diagram to the I 
coefficient is given by (for both % — V, A) 

p{d) _ TU )2 [ d D l E^jv ■ k)k + + E 7 (/ + fc + - k ± ■ l ± ) + \L(k ± ■ h - l+k+) 

H n) J {2<k) d (-v I + ie){{l + k - Pl ) 2 - m 2 + ie){l 2 + 21 ■ k + ie){l 2 + ie) 1 ' 

Here the contribution of the 0(l 2 ) terms is of order 1/-E 7 and thus subleading. The first 
two terms, ~ 1° and ~ I 1 , can be computed to leading order in E 1 by expanding the large 
denominator as (I + k — p 7 ) 2 — m 2 ~ — 2E 1 (l + + k + ). The numerator of the first two terms 
can be arranged as the sum of two terms, one of which just cancels the denominator l + + k + , 
plus a remainder 

E,,2{v ■ k)k+ + E,,(l + k + - k± ■ l±) = £ 7 U{v ■ k)(l+ + k+) + ~ k ± • l S) • ( 57 ) 

The first term has exactly the structure of the scalar integral appearing in the correction to 
/b- We obtain for the total contribution of the box diagram to leading order in E 1 as 

£W = iUn) 2 1 2(v ■ k)J m (v ■ k) + / Srni ; w , /+ ^+" ~ fc ^ r^ 7775 ~\ (5* 

v ; 1 v ' K J J {2ix) D (-v ■ I + ie){l + + k + )(l 2 + 21 ■ k + ie){l 2 + ie) } y 

with Jm(ko) the IR singular integral defined and computed in the Appendix (see Eq. flAllp ). 
The second integral in ([58]) is IR finite and can be easily computed by combining the last 
two denominators with Feynman parameters. This gives 

in d + M _ m 



(2vr) 4 (-vl + ie)(l+ + k + )(l 2 + 21 ■ k + te)(l 2 + ie) 
1 r dH (l + — xk + ,l± — xk±) 

' r (2tt) 4 (-vl + xk + ie)(l + + fc+(l - x)){l 2 - m 2 x 2 + ie) 
"dxTdL (! + -xk + ,-xk ± ) 



2 



8tt 2 Jo Jo + [l + + Jfe+(1 - x)} [l\ - 2xl + k + m 2 x 2 - ie] 

where we performed the integration over the light-cone coordinates l-,l±. Inserting these 
results into the second integral in Eq. Q5BD one obtains 



1 {k±)2 [1 dxHdl + - — — 9 l+ - , — (60) 



^7r 2 k + Jo Jo [/+ + k + (l — x)][l\ — 2xl + k + m 2 x 2 — ie 
log 2 — - + 2ni log ; 



(4tt) 2 V 2k Q to 2k J ' 

Finally, we obtain the following result for &^ to leading twist 



= f$ = z(4tt) 2 2k J IR (k ) + log 2 |±- + 2<Ki log ^- . (61) 
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Diagram 


Contributions to £i(E 7 ) 


2(a) 
2(b) 
2(c) 
2(d) 

\{zr°-i) 

1 {Z HQET_ 1) 


N uv log ^± i og ^ + i og ^ 
N uv log 2fe +2 log 2 ^ 2 fc+ 1 
N uv + log 2S 7 fc + ! 

2(4vr) 2 2k J IR (k ) + log 2 ^ + 2ttz log ^ 

^ N UV + 3 log ^_ N IR 
N UV _ N IR 


Total 


3 N UV log 2^+ + l log 4 

1 2 2^ . i 2E- 1 . i 2E-,k+ 4^2 

log fc ; +io g fc ; +2 log ^ + 3 

+A /R (A; )+log 2 ^ + 27r2log^ 



Table 1. One- loop contributions to the form factor from individual diagrams. 
The IR singular contribution A IR (k ) = i(4n) 2 2k J IR (k ) - 2(N I £ R + log^r) is 
identical to the one appearing in the one-loop correction to /# fl23D and can be 
absorbed into the B meson light-cone wave function as explained in Sec. II. 

We are now in a position to write down the complete one-loop correction to the form 
factors for B — > ip£j. The individual contributions from the diagrams of Fig. 2 and their 
total result are presented in Table 1. There are a few remarks which can be made about 
these results. 



The box diagram (Fig. 2(d)) contains an IR divergent term (|61"D which depends on 



k± through the quantity v ■ k. Note that this is different from the case of the pion 



form factor, which is IR finite and contains only mass singularities. However, 
the IR singular term can be seen to be precisely identical to the one appearing in the 
one-loop correction Eq. (|22|) to the decay constant fg. As explained in Sec. II, it can 
be absorbed into the B meson light-cone wave function, leaving a IR-finite Wilson 
coefficient depending only on the light-cone momentum component k + . 



• The dependence on the MS hybrid scale /ih cancels, as it should, between the £ coeffi- 
cient and the corresponding factor in (|42]). We will choose for this scale fi = 2E 1 , with 
which the first factor accounts explicitly for the large logs (^-log(^-)) n in leading 
logarithmic approximation. 

• The equality of the leading twist form factors for different currents noted at tree level 
fv(E^) = /a^E^) persists to one-loop order. In view of the symmetry arguments 
justifying this equality at tree level, it is tempting to conjecture that this is a general 
result for the leading twist form factors, valid to all orders in the strong coupling. 

With these remarks, the leading twist result for the form factors in B — > Ivf) decays can 
be written asQ 



4 The modified wave function ?/>(£;+, kj_) in this expression contains only the exponentiated IR 
singularity. 
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where the hard scattering kernel Tn{k + , k±) is given to one- loop order by 

T H (k + , k ± ) = — ( 1 + sV ,_ 7; <| - log 2 + - log + - log — ^± 



£;+ \ 47T I A; + 2 fc + 2 m 

47T 




+ lor [l + j ^|-27rilog(l + 7 i^| . (03) 

Note that this result for the form factors is sensitive to the dependence of the wave function 
on the transverse momenta, through the last two terms. After integration over kj_, these 
terms will give a finite correction to the light-cone wave function. The last term in (|63|) will 
give the form factors also a complex phase. These features are in contrast to the pion form 
factor case, where transverse momentum effects are absent to leading twist. 



IV. MASS-SINGULAR LOGARITHMS AND SUDAKOV EFFECTS 

The expression for the form factors ( |62j) contains mass-singular logarithms 
log(2i? 7 A; + /m 2 ) as well as Sudakov double logs log {2E 1 /k + ) which must be resummed to 
all orders. In this section we will discuss these issues in turn. 



A. Resummation of collinear singularities 

The hard scattering amplitude Tn{k+) contains collinear logs log(2i? 7 A; + /m 2 ), which 
arise only from the diagram 2(b) and the wave function renormalization constant. In the 
former, these logs are produced by integration over the transverse momenta in the region 
m < 5. < 2E^k + . The propagator of the struck quark can be written in this region as 
(k— p-y + l) 2 — —2Ery(k + + l + ) + - ■ •. Keeping only the leading terms in l/E^, the contribution 
of this diagram (plus the tree contribution) is proportional to 



p(o) + r (6) _ 



(64) 



(k, 



2k, 



log 



2E^k + 
m 2 



J dk + di + i>(k + )ls(k + -i + )+)c(k + ,i + ) J m2 -j±\t£\i + ), 

where T^\l + ) = is the tree-level hard scattering amplitude. The kernel JC(k+,l+) is 
given by 



)C(k + ,L 



4:71 



' 2i + e(k + -i + ) \ i " 

k + (k + -l + ) )+2 6 ^-^ 



(65) 
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where the term proportional to S(k + — l + ) comes from the wave function renormalization. 
The +-distribution is defined as usual by 

/(*+, = f{k + , 1+) - 5(k + - l + ) J dr + f(k + , r+ ) . (66) 

Integrating over l + gives the explicit one- loop result for the mass-singular logs (|53|) . 
However, writing the result in this form helps us to resum these logs to all orders. To do 
this, one cuts the integral over the transverse loop momentum in ( |64"D to a certain cut-off 
A. This will be chosen identical to the one introduced in the normalization condition (P7|). 
The logarithm resulting from integration over the range m 2 < f]_ < A 2 is then absorbed into 
the wave function ip{k + , A) by defining 

4>(l + , A 2 ) = J dk + ^(k + ) (^5(k + - 1+) + K(k+, 1+) log ^) . (67) 
Expressed in terms of the wave function ip(k + ,A), the form factor is written as 



Taking the logarithmic derivative of (57 ) with respect to A 2 gives the integral equation 
satisfied by if)(l+, A 2 ) 

A 2 — A 2 ) = / dk + JC(k + , l + )i>(k + , A 2 ) . (69) 



dA 2 " v ' Jo 

This is the analog of the Brodsky-Lepage evolution equation which governs the evolution of 
the light-cone wave function of a heavy meson with the factorization scale. The moments 
of the wave function are renormalized multiplicatively with the anomalous dimensions 

where (k 1 !) = J dkj r (k + ) n il)(k + , A). The th moment of the wave function evolves with the 
same anomalous dimension 3a s (7F/(87r) as previously derived in Sec. II (see Eq. ( p7|) ). 



B. Sudakov resummation 

The radiative correction to the form factors fi(E^) contains double logarithms of the large 
ratio log 2 The explicit calculation of the preceding section shows that such logarithms 
arise (in the Feynman gauge) from one loop correction to the vertex b — ► Wq of the weak 
decay of the h— quark into a light quark with momentum Q = p 7 — k. It is easy to see that 
in the rest frame of the B— meson in the kinematical region E 1 ^> Aqcd the light quark 
moves close to the "+" light-cone direction along the photon momentum with the energy 
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(Q ■ v) ~ £ 7 and small virtuality Q 2 2(p 7 • k) such that Q 2 /{2Q ■ v) 2 = O (A/£ 7 ). 

It is the ratio of the scales Q 2 /(2Q-v) 2 <C 1 that enters as an argument into Sudakov 
double logs. The appearance of large negative corrections is related to enhancement of the 
contribution of soft virtual gluons propagating collinear to the produced light quark, close to 
the direction of photon momentum. In contrast with the inclusive distributions where it is 
canceled against the contribution of real soft gluon emissions, virtual soft gluon contribution 
survives for an exclusive distribution like the one under consideration due to the absence of 
real soft gluons in the final states. 

Let us consider the one-loop Sudakov correction to the weak decay vertex 

F=l-^C F S, S^-^f-X^-^f (71) 

with 2(Q ■ v) = 2Ery and Q 2 = —2E^k + in the rest frame of the B-meson. The Sudakov form 
factor S is given by the following one-loop Feynman integral 

* = (4Q_ - 2Z_K 

J (2vr) 2 (-v ■ I + ie){l 2 + ie)((Q - I) 2 + ie) 1 ' 

with Q_ = 2E y , Q + = —k + , Q 2 = Q+Q- and v + = 1. Calculating S we write the integration 
measure as d A l = ^dl + dl-d 2 l± and perform l + — integration by taking residues at the poles 
corresponding to three different propagators. One finds that the integral is different from 
zero provided that /_ belongs to one of the regions, < /_ < Q- and /_ > Q_. One checks 
that in the second case the gluon has large components of the momenta and its contribution 
is associated with short distance (hard) subprocess. In the first case, the /+— integral is given 
by the residue at l + = l\/ (2Z_) — ie which effectively amounts to putting the virtual gluon 
on-shell, I 2 = 0. Then, introducing the scaling variable x = 1-/Q- one finds 

S 4/-f<7 ^. (73) 

The denominators effectively set the limits on the integration ranges, such that the leading 
doubly logarithmic correction arises from the region (in the rest frame of the B meson) 

k+ < I- < 2£ 7 , k + l_ <l\<l 2 _. (74) 

In this way, one calculates the one-loop correction to the Sudakov form factor as 



k+ l- Jk+i- lj_ \ 2E, 

2k+Ei dl ± ,„ l ± , [ {2Ei? <ti± ,„ ( 2E -y) 2 r {2E ^ ? dl ± 

7 l± 



7 > 

r> ln 7T + /' ' ^ln^?^- / ^ (75) 



The reason why we represented the one-loop correction in this particular form is that it ad- 
mits generalization to higher orders in the coupling constant that effectively resums Sudakov 
logarithms. Each term in the r.h.s. of (|75|) comes from different part of gluon phase space 
and has the following interpretation. The last term describes collinear emission of on-shell 
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energetic gluon, L = 0{E 1 ), l + <C Z_ and l\ = /+/_ = 0(k + E y ) (collinear region), and gives 
rise to a single collinear log. The first two terms correspond to soft gluon emission on two 



different infrared scales, l + ~ L ~ l± = 0{k + ) (soft region) and l + ~ /_ ~ l± = 0{^k + E 1 ) 
(infrared region), and produce double logarithmic contributions. Examining higher order 
corrections to the b — > Wq vertex one can show that the same regions of gluon momenta 
provide the dominant contribution to the Sudakov form factors. Moreover, since the emis- 
sion of collinear and soft gluons occurs on different time scale their contribution factorizes 



out as |T!J 



F = F H (Q 2 + , fi 2 )Fj(Q 2 + , Q + Q-,fx 2 )F s (Q 2 + , Q + Q_, fx 2 )F m (Q + Q_, Q 2 _,fi 2 ) (76) 

with Q_ = 2E 1 and Q + = —k + . Here the hard subprocess Fjj takes into account short 
distance corrections to the weak decay vertex, ~ Q + , while Fj, F s and Fir denote contri- 
butions of collinear and soft gluons on different momentum scales. The parameter [i entering 
([76]) plays the role of the factorization scale. The subprocesses Fs and Fir admit an operator 
definition as expectation values of the Wilson lines originating as eikonal phase of quarks 
interacting with soft gluons. Using this interpretation one can show that the fi— evolution of 
Fs and Fir subprocesses is in one-to-one correspondence with the renormalization proper- 
ties of Wilson lines. In this way, using evolution equations for different subprocesses and the 
[i— independence of F, one finds that the Sudakov form factor obeys the following evolution 
equation 

dlnF „, , „ , 1 /-(2^) 2 dl 2 . 



d\nE y 



i r v 2E i) dr 

T(a s (Q 2 + )) + T (a s (Q 2 _)) - - / ^Wa.(£)) . (77) 

I J2E^k+ It 



This evolution equation involves three functions of the coupling constant that appear as 
anomalous dimensions in the evolution equations for different subprocesses. Two of them 
r cusp and To are related to renormalization of (light-like) Wilson loops while T is related to 
the UV renormalization of the weak decay vertex 

r cusp = -C F + 0(a 2 ) , T = + O(a 2 ), T = —Cf + 0(a 2 ) . (78) 

7T 7T 



Neglecting the r — term one can write the solution to the evolution equation Cf77Q as 
' Q2 dl\, 11 ^ , , f^dP, Q* , .... rQldP 



41nF= / -^ln-^r cusp (a,(/i))+ / -± In ^±r cusp (a s (/i)) - / -±T(a s (l 2 ± )) 

JQj_ L ± JQ 2 (j_ l ± JQ 2 L ± 

(79) 



Comparing (|79|) with the one- loop expression (|7"5|) we conclude that, first, Sudakov loga- 
rithms exponentiate f22| and, second, the exponent of the Sudakov form factor is formally 



given by the one-loop expression in which the "bare" QCD coupling constant is replaced 
by an anomalous dimension with a particular choice of the normalization scale given by 
gluon transverse momentum l\. The perturbative expansion ( |79|) is valid provided that the 
integration over l\ does not go below the Landau singularities of the coupling constant. 
This means that the resummed expression ([79]) is valid provided that k 2 + 3> Aqcd- in the 



practical application discussed below, the Sudakov form factor for k + below the singularity 
at k + = Aqcd will be frozen at its value just above this point. 
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Expanding the QCD running coupling constant ce s (l\) in powers of a s {E 1 ) and per- 
forming the integration in (|79|) one can expand the exponent In F into a series of the form 



[a s L 



2\n 



Of' 



with n = 1 , 2 



a. 



a s (i? 7 ) and L = ln(E 1 /k + ) to which we 



shall refer as leading-order (LO), next-to-leading order (NLO), . . . corrections. In particular, 
to the LO approximation it proves enough to keep only the first two terms in (|79|). Using 
the one-loop running of the strong coupling 



47T 



A)hi|l 



(80) 



with j3 = 11 — |nj and replacing r cusp by its one-loop expression, one gets 



S 



2C F \ 2E 1 k + 
i,o = -^|-log^log 



1 , 2E^k+ 

— log ■ 

2 8 A 2 



2Ery 2E*y /C-(- 

+ log log log — + log — log log — 



(81) 



which agrees with the result in |20|. The overall effect of the Sudakov form factor is to depress 
the form factors at large values of E^. One can systematically improve the accuracy of ( |79| ) 
by taking into account NLO terms. To this end one should include two-loop corrections to 
the coupling constant and r cusp as well as one-loop correction to the anomalous dimensions 
T defined in (fl£\). 



V. APPLICATION 

The decay rate for B — > 7^£ + differential in the lepton and photon energy is 
_^T_ = *GWutfrn% + f 2 {Ei)][ _ 2xy + 2xy 2 + x _ 2x * y + ^ 

-2f A (E 7 )f v (E 7 )x(l -x)(l+x- 2y)} . 



(82) 



We denoted here x — 1 — 2E 1 /rrtB and y = 2E e /rriB, in terms of which the available phase 
space is described as x = (0, 1) and y = (x, 1). An integration over all possible values of the 
electron energy y gives for the rate as function of the photon energy 

dL _ aG%\V^m%_ [fl{E ^ + fv{Ej)]x{1 _ x)S _ (g3) 



dE y 12(2vr) 

Our results for the form factors fy^iE^) can be therefore turned into a prediction for the 
shape of the photon spectrum in this decay. To leading twist, the 1/E^ dependence of these 
form factors yields a symmetrical photon spectrum dT/dx oc x(l — x). 

Neglecting radiative corrections, the form factors parametrizing the B + — > £ + i"y decay 
are given by 

a^>-^<<w-£> +0 (£), m 

where we included also the leading A/m^ correction computed in (p5|). 
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Extrapolating the tree-level form factors (|84D over the entire phase space gives for the 
integrated decay rate 



This result is identical to the one obtained in |JJ from a quark model calculation of the 
annihilation graph, with the identification R — > l/m u (the inverse constituent quark mass). 
In fact the appearance of the inverse constituent quark mass is a common aspect of quark 
model calculations of long distance effects produced by weak annihilation topologies with 
emission of one photon or gluon ||25|1 . Such contributions have been investigated in many 



processes such as B — > p7 |26| and B — > D*"f 



Our QCD-based derivation gives such computations a precise meaning by replacing the 
ambiguous notion of constituent quark mass with a well-defined integral over the light-cone 
B meson wave function. Besides specifying the limits of validity of this result, such an 
approach allows one to compute also strong interactions corrections to it in a systematic 
way. 

It is possible to derive a model-independent lower limit on the magnitude of the R 
parameter, under the assumption that the light-cone wave function is everywhere positive, 
which is reasonable for the ground state B meson. This bound reads 

**8 = sc- (86) 

and can be proved with the help of the inequality 

-?- + bk + >2Vb, k + >0. (87) 
k+ 

Here b is an arbitrary real positive number. Multiplying with ip(k + ) and integrating over k + 
gives the inequality R > — |A6, where we used the normalization condition ([H]). This 
is most restrictive provided that one chooses yb = jr, which gives the result (|36|). 

It is interesting to note that a hadronic parameter related to R appears also in the 
description of the nonfactorizable corrections to nonleptonic B — > tttt decays (called 



there 1/A^). Our results suggest therefore a method for extracting this parameter in a 
model-independent way from data on B — > jeu decays. 

To eliminate the dependence on fs and V u b, we will present our results for the photon 
spectrum by normalizing it to the pure leptonic decay rate for B + — > fiu, which is given by 

r^ + ^ H = ^|« / l(^) 2 ( 1 -4). (88) 

8n \m B J \ m B J 

For illustrative purposes we will adopt in the following numerical estimates a two- 
parameter Ansatz for the heavy meson light-cone wave function inspired by the oscillator 
model of |24] 



^(k + ) = ATk+exp ( -^#+ -a) 2 ) . (89) 
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We will vary the width parameter oj in the range to = 0.1 — 0.3 GeV. The parameters TV and 
a will be determined from the normalization conditions discussed in Sec. II. For a given value 
of A, these normalization conditions set an upper bound on the width parameter u, given 
by uj max = -^t^ A (corresponding to a = 0). The latter will be taken between A = 0.3 GeV 



and 0.4 GeV. The resulting numerical value of the constant R together with the parameter 
a are given in Table 2 for several choices of A and u. 



u (GeV) 


A = 0.3 GeV 


A = 0.35 GeV 


A = 0.4 GeV 


0.1 


0.2 


0.3 


0.1 


0.2 


0.3 


0.1 


0.2 


0.3 


a (GeV) 


0.37 


0.27 


0.05 


0.44 


0.36 


0.19 


0.51 


0.44 


0.30 


R (GeV" 1 ) 


2.70 


3.29 


3.87 


2.28 


2.65 


3.09 


1.96 


2.23 


2.59 



Table 2. Light-cone wave function parameters a and R corresponding to several 
values of the binding energy A and the width parameter u. 



Taking f B = 175 MeV and \V ub \ 
mode a branching ratio 



3.25 x 10 3 gives for the muonic B + decay 



B(B + -> fiu) = 2.3 x 10" 



(90) 



For a typical range of values R = 2 — 3 GeV 1 (see Table 2), the tree- level integrated rate 
(jS5|) predicts a ratio 



B{B+ -> 7 e+z/) 
B(B+ -> fiu) 



2.02R 2 



(91) 



which implies branching ratios of about (2 — 5) x 10~ 6 for the B + radiative leptonic mode, 
in agreement with the general estimates of ||. 

A similar analysis can be made for the radiative leptonic decay D + — > 7e + z/, for which 
one obtains the ratio of branching ratios 



B{D+ fiu) 



■ (),,R- — 



0.09-0.16 



(92) 



Note that the charmed quark contribution can be appreciable, and can account for up to 
50% of the light quark contribution. Neglecting SU(3) breaking effects and small kinematical 
corrections, the denominator can be related to the muonic branching ratio for decay 
which has been measured 



B{D A 



Vcs 



-B(D 



(0.68 ±0.37) x 10" 



(93) 



We used here the CLEO result B(D S -> p+v) = (6.2 ± 3.1) x 10" 3 0. This predicts an 
absolute branching ratio for the radiative D + decay of 



B{D- 



7 e+z/) = (0.82 ±0.65) x 10" 



(94) 



Somewhat larger absolute values are obtained for the radiative decay width, which is 
enhanced by the larger CKM matrix element V cs . Neglecting SU(3) breaking in the hadronic 
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parameter R one finds for this case B(D+ — > 7e + z/) = (0.9 ±0.8) x 10 -3 , again in agreement 
with the estimates of ||. 

While useful as an order of magnitude estimate, we stress that the relation (|5|) and the 
numerical results obtained with its help are not rigorous predictions of QCD in any well- 
defined limit. The reason for this is that the prediction (Q) for the form factors fv,A{Ey) 
receives uncontrollable corrections of order A 2 /E 2 as soon as the photon energy E 1 does 
not lie within the region of applicability of our analysis Aqcd "C E t A similar statement 
can be made about the corresponding predictions for the charged lepton energy spectrum, 
which requires knowledge of the form factors over the entire range of E T 

In order to avoid these problems, we will restrict our considerations to quantities defined 
with a sufficiently high lower cut on E 1 . When radiative corrections are taken into account, 



the hadronic matrix element R in (|84| ) acquires a logarithmic dependence on E 1 given by 



W ' f^)^7^ tt %l( 1 + ^ W ) +0 f*i) . (95) 
v 11 \a s {2E 1 )J Jo + k+ V 4vr v 7 V \E^ m b J K ' 

We show in Fig. 3 the results obtained for the form factors and in Fig. 4 for the photon energy 
spectrum using the tree-level form factors and including the one-loop correction computed in 
Section III. This correction decreases the rate, at least in the region of validity of our results. 
This effect is mostly due to the double log in the one-loop hard scattering amplitude; the 
leading-log factor in (^) makes a positive contribution. This illustrates the importance of 
the double logarithms log 2 (^ 1 ), which have to be resummed to all orders. The third curve 
in Figs. 3 and 4 shows the spectrum obtained by resumming the Sudakov logarithms to all 
orders, as explained in Sec. IV. 

While the functional form of the hadronic matrix element R{E 1 ) depends on the detailed 
form of the (unknown) B meson light-cone wave function, it is important to note that it 
is independent of the heavy quark mass m b (up to calculable logarithmic corrections). One 
would like to eliminate it by taking ratios of the photon spectra in B and D radiative leptonic 
decays. However, the large value of the l/m c correction in the latter case would introduce 
large corrections to such a ratio, which shows that some knowledge of R{E 1 ) is necessary. 

With this view in mind, we propose in the following a two-step procedure for determining 
the magnitude of the CKM matrix element \V ub \. In the first step, the hadronic function 
R( b '(E y ) is determined in a region ^> Aqcd from the normalized photon spectrum in B + 
decays 



1 dT(B + -> 7e+z/) _ am B / R (b)( E ) Qb\ ( m s\ x B (l 



T e (B + — > dE 1 37T V m b J \m M I 1_ 



2 



III 
III- 



(96) 



D 



with xb = 1 — (2-E , 7 )/m^. We used on the RHS the leading twist result for the B + form 
factors; the 1/m.f, correction is very small and will be neglected. The superscript on R^ b \E^) 
labels the heavy quark flavor. 

In the second step, one takes the ratio of photon spectra in B and D decays, which is 
given by 



V ub A ( Q u R {b \E, 



V r 



cd 



(97) 



^R^iE^-QJm 
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V, 



nb 



cd 



\ QdCR {b) ( E -f) - Qc/mJ \a s (m c )l \m D ) x D 



where R^ b \E 7 ) is known from (|96|). We used here the logarithmic dependence on the heavy 
quark mass © for the R^ Q \E 1 ) coefficients 

(/ \ \ -2//3 
aJm b ) I 



and the large mass scaling law fL4]|[5| for the pseudoscalar decay constants 



f B _ f a s (m b ) \- 2,Po /W 
f D \a s (m c ) V m B 



(99) 



The result (|97|) can be used to determine the CKM matrix element \V ub \. 

The leading corrections to this determination come from higher-twist effects of order 
0(A 2 /E 2 ) in the D + meson radiative leptonic form factors. Their magnitude can be esti- 
mated by comparing the normalized photon spectra ([56|) in B and D decays. Although for 
the B case these corrections are expected to be well under control over a reasonably wide 
range of values for E y , it is questionable whether in the D case such an large energy region 
Aqcd E 1 exists at all. Since the maximum photon energy accessible in D decays is only 
about 0.93 GeV, the higher twist effects can be expected to contribute no less than 10% to 
the D meson form factors. 

A similar determination of \V ub \ can be performed using instead of D + , the more ac- 
cessible Df meson radiative leptonic decays. However, this would introduce an additional 
uncertainty on the theoretical side through SU(3) breaking effects. 



VI. CONCLUSIONS 

We studied in this paper the form factors for the radiative leptonic decay of a heavy 
meson (e.g. B + — > r ye + u) in an expansion in powers of the inverse photon energy 1/E^. 
To leading order 0(1/ E y ) these form factors are given by a convolution of the light-cone B 
meson wave function ip(k + ) with an infrared- finite hard scattering kernel T^(fc + , k±) (j63|). 

Physically, this problem is very similar to the Tf-f-f* pion form factor F nl (Q 2 ) studied in 
IjTiyiBfl , where a similar factorization can be established for the leading twist contribution 
of 0(1/ Q 2 ). However, there are some important differences, the most striking of which 
concerns the dependence on the transverse momentum to leading twist revealed in the form 
of the hard scattering amplitude Tn(k + , k±). Such a dependence is absent in the case 
of the pion form factor, and its appearance can be traced to the presence of the second 
dimensional parameter k + (the light-cone projection of the light quark momentum in the B 
meson) in addition to the large scale E 1 . On the practical side, this implies a certain loss of 
predictive power: while the logarithmic dependence on E^ is well-determined, the constant 
term depends on the precise form of the full 3-dimensional light-cone B wave function. 

A second important complication compared to the F^(Q 2 ) case consists in the appear- 
ance of Sudakov double logarithms, which have to be resummed to all orders. This feature 
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has been noted previously in the context of the B — > 7r(p) semileptonic form factors of a 
heavy hadron in [2^.21], where these Sudakov effects have been resummed (up to next-to- 
leading order). Numerically, their effect is most important near the upper end of the photon 
energy spectrum. 

An interesting qualitative result of our analysis is the equality of form factors of different 
currents fy(E 7 ) = /^(£? 7 ) at leading twist. While this equality was established by an explicit 
one-loop calculation, it is probably a general result, true to all orders in the strong coupling. 
In a perturbative QCD language, the reason for this equality roots in the dominance of the 
momentum integration regions where the propagator of the struck quark (see Figs. 2) can be 
approximated with a light-like eikonal line. This relation can be formalized by going over to 
an effective theory jl7f where the couplings of gluons to this line possess a higher symmetry. 

A similar approach has been taken in |TJ| to derive relations among semileptonic decay 
form factors of a heavy hadron. However, in the latter case the hard one-gluon exchange 
mechanism can be shown to introduce corrections to these relations, already at leading 
twist. This is different from our case where these relations appear to be preserved (at least 
at one- loop order) under inclusion of the hard gluon exchange. 

Finally, our formalism can be used to put previous quark model estimates of radiative 
leptonic decays on a more firm theoretical basis, by giving a precise definition of the 
light quark constituent mass. Our approach is likely to give a reliable description of the 
form factors in the large E 1 region, up to controllable corrections of order A 2 /E 2 . This 
complements an alternative approach presented in |J which is best suited to the low-E^ 
region, where the heavy hadron chiral perturbation theory is expected to be applicable. 

Using as input parameter the binding energy of a B meson A, we gave several estimates 
for the branching ratios of these modes. As a by-product, we presented also a method for 
extracting the CKM matrix element \V u b\ by comparing photon energy spectra in radiative 
leptonic B and D decays. 
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APPENDIX A: SCALAR INTEGRALS 



We present here a few details relevant for the computation of the radiative corrections. 
The scalar integral appearing in the heavy-light vertex correction is computed by first 
combining the two massless propagators with the help of a Feynman parameter 



d 4 Z 1 



(2tt) 4 (-vl + ie)((l + q - k) 2 + ie)(l 2 + ie) 
•i , r d 4 / 1 



(2tt) 4 (-u • Z + ie)((Z + x(g - A;)) 2 - s + ie) 



2 



/ d/ + / d;r 79 ^1 , ( Al ) 

8tt 2 7o + Jo l\- 2xq l + + s v 1 

with s = — x(l — x)(q — A;) 2 . After shifting the loop momentum Z — > Z — x(g — fc), one 
integrates over the light-cone component Z_ using the Cauchy theorem, and subsequently 
over the transverse momentum Zj_. 

The integral with one power of Z Q in the numerator can be reduced to a two-point function 
plus a UV finite integral by first combining the massless denominators with a Feynman 
parameter as above. This gives for the numerator 

L -> L - x(q - k) a = (v ■ l)v a + (Z_l) q - x(q - k) a (A2) 
= -[-v • Z + xv ■ (q - k)]v a - x[(q - k) a - v ■ (q - k)] + (l±) a . 

The first term cancels the heavy quark propagator in the denominator, and the (l±) a term 
vanishes after integration over I. One obtains in this way 

/(«) = f — hi. (AD 

a ~ J (2tt) 4 (-v ■ I + ie)((l + q - k) 2 + te){P + ie) 1 J 

= - Ua £« - [(q -k) a -v(q- k)](xJ^) (A4) 



with 



V(«) = / I 

J (2vr)™[(Z + g- k) 2 + te}{l 2 + ie) 1 ; 



and 



(xjW) = ^ / dZ+ / drr-, =- (A6) 

These integrals can be evaluated exactly with the following results 



(47r) 2 (E 7 -A) v ^\ Z V 2S 7 Mv^+l)/ V 2^(^+1) 

«i - ,)^) = » - M">> = i ( i + i) j + {° log i^f^ 
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+ (1 + 

— 6 log 
+ 



n 8(E,-A)Vg 



2£ 7 fc+ 
2v^ 



; & 2£ 7 A; + + r 2 



2£ 7 £; 



2£ 7 £; + 



log 



2.ELfc + + T X 



1 + (E 7 -A) 2 (l + t, 2E 1 k + + r 2 

{2E 1 k + ) 2 2E 1 k + + n 



log- 



(£ 7 - A) 2 (2£ 7 A; + - r a )£ b 2£ 7 £; + (l + VVf) 

{2E 1 k + f 2E 7 k+ + r 2 \ 

(£ 7 - A) 2 (2£ 7 Ac + - r 2 )£ ° g 2E 7 A; + (1 + 1/^) J 

with £ = 1 + 2E 1 k + /(E 1 - A) 2 and r lj2 = 4(£ 7 - A) 2 (±l + y^) ± 2£ 7 Ac+. 

The vertex correction to the photon coupling to the light quark is parametrized in terms 
of the integrals 
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iVf^ + 3-log 



2p 7 • fc' 



When computing the one-loop correction to fs and the box diagram, one encounters the 
IR singular integral 



Jir(v •k)=fi 



2t 



d U l 



1 



(2tt) d (-vl + ie)(l 2 + 21 ■ k + ie)(l 2 + is) 
(47r// 2 rr 2 ) e 



(A10) 



(4tt) 



-r(l + e) / dx 
e Jo 



[(xk - v) 2 - ze] 1+e 



The IR singularity has been regulated with dimensional regularization in D = 4 — 2e dimen- 
sions. The integral over x can be computed explicitly with the result (with Q = v ■ k) 



Jir{Q) 



+ Li 2 



(4vr) 2 2^Q 2 - m 2 
' 2^ 



{N* R + log 



/i- 



log 



Q - VQ' 2 - m 2 



Q + VQ 2 ' 



2ixi 



(All) 



- Li 2 



-2^ 



Q-VQ 7 ' 
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~6~ 



2i\i log 



4(Q 2 - m 2 ) 



77T 



with N£ R = 1/e — 7s + log(47r). In the limit Q S> m this agrees with the expression given 
in the Appendix C of 
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FIGURES 




FIG. 2. One-loop corrections to the radiative leptonic decay B + — > W7. The curly line 
represents a gluon. The quark wave function renormalization corrections are not shown. 
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FIG. 3. Typical leading twist form factors fi(E^) [i = V,A) for B — > ^ev decays. The 
continuous line shows the tree-level result, the dotted line includes one-loop corrections to the hard 
scattering amplitude, and the dashed line includes the resummed Sudakov form factor truncated 
with a cut-off at (k + ) m i n = Aqcd- We use a s (rrib) = 0.3 and a = 0.36 GeV, oj = 0.2 GeV, 
corresponding to A = 0.35 GeV. 




E (GeV) 

FIG. 4. The photon spectrum in B + — ► je + i / normalized to the pure muonic leptonic decay 
rate. The continuous line represents the tree-level result, assuming R = 3 GeV -1 . The dotted 
line includes the effects of the one- loop strong correction (only the logarithms) with a s (m,b) = 0.3, 
and the dashed line includes the resummed Sudakov logs. The same parameters are used for the 
light-cone B wave function as in Fig. 3. 
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